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Abstract. Various properties of the Misner- Sharp spherically symmetric gravitational 
energy E are established or reviewed. In the Newtonian limit of a perfect fluid, E yields 
the Newtonian mass to leading order and the Newtonian kinetic and potential energy 
to the next order. For test particles, the corresponding Hajfcek energy is conserved and 
has the behaviour appropriate to energy in the Newtonian and special-relativistic limits. 
In the small-sphere limit, the leading term in E is the product of volume and the energy 
density of the matter. In vacuo, E reduces to the Schwarzschild energy. At null and spatial 
infinity, E reduces to the Bondi-Sachs and Arnowitt-Deser-Misner energies respectively. 
The conserved Kodama current has charge E. A sphere is trapped if E > hr, marginal if 
E = \r and untrapped if E < |r, where r is the areal radius. A central singularity is spatial 
and trapped if E > 0, and temporal and untrapped if E < 0. On an untrapped sphere, E 
is non-decreasing in any outgoing spatial or null direction, assuming the dominant energy 
condition. It follows that E > on an untrapped spatial hypersurface with regular centre, 
and E > ^tq on an untrapped spatial hypersurface bounded at the inward end by a 
marginal sphere of radius tq. All these inequalities extend to the asymptotic energies, 
recovering the Bondi-Sachs energy loss and the positivity of the asymptotic energies, as 
well as proving the conjectured Penrose inequality for black or white holes. Implications 
for the cosmic censorship hypothesis and for general definitions of gravitational energy are 
discussed. 
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I. Introduction 

A massive source produces a gravitational field which has energy. In Relativity theory, the 
equivalence of mass and energy means that it is only the combined energy which may be 
measured at a distance. Moreover, the non-linearity of the gravitational field means that 
the material (or passive) mass and its gravitational and kinetic energy combine in a non- 
linear, non-local way to produce the effective (or active) energy. In spherical symmetry, in 
vacuo, this effective energy is just the Schwarzschild energy. In general, there is no agreed 
definition of this energy, except at infinity in an asymptotically flat space-time, where 
one has the Arnowitt-Deser-Misner [1] energy Eadm at spatial infinity and the Bondi- 
Sachs [2-3] energy Ebs a t null infinity. One would therefore like a definition of energy 
which reduces to these asymptotic energies appropriately. Also, given the above physical 
motivation, one would like the energy to yield the Newtonian mass in the Newtonian 
limit, with the highest-order correction yielding the Newtonian energy. Similarly, one 
would like the energy to yield the correct energies of test particles in the Newtonian and 
special-relativistic limits. Remarkably, such an energy does exist in spherical symmetry: 
the Misner-Sharp [4] energy E. Moreover, E is intimately related to the characteristic 
strong-field gravitational phenomena, namely black and white holes and singularities. 

This article lists various key properties of this energy. The main new results, apart 
from the Newtonian, test-particle and special-relativistic limits, are a monotonicity prop- 
erty of E which leads to a positivity property and a lower bound for E in terms of the 
area of a black or white hole, the so-called isoperimetric inequality [5]. Here, black and 
white holes are defined by marginal surfaces of certain types, as explained later, or in [6] or 
[7] . These properties of E extend to the asymptotic energies, in particular establishing the 
isoperimetric inequality for Eadm, as conjectured by Penrose [5], and for Ebs- Known 
results are also reviewed, partly because some are prerequisites for the newer results, and 
partly because the existing literature on the subject is somewhat dispersed, so that it is 
not always appreciated that E enjoys quite so many desirable properties. 

In Section II, E is defined geometrically and shown to have various purely geomet- 
rical properties related to trapped and marginal surfaces, central singularities and the 
asymptotic energies. In Section III, various dynamical properties are derived assuming the 
dominant energy condition, including the monotonicity, positivity and area-bound proper- 
ties. A discussion of implications for the cosmic censorship hypothesis is also given there. 
In Section IV, the geometry is decomposed with respect to spatial hypersurfaces and the 
behaviour of E in the Newtonian limit is found. In Section V, the energy of test particles 
as determined by E is discussed, including the Newtonian and special-relativistic limits. 
The Conclusion discusses the implications for more general suggestions for gravitational 
energy. The Appendices concern (A) Newtonian theory, (B) the identity of E as a charge 
associated with a conserved current, and (C) energy-momentum. 

Since the whole article is concerned with spherical symmetry, this case will be assumed 
implicitly in the Propositions, without repeated qualification. Similarly, all geometrical 
objects mentioned will be assumed to respect the spherical symmetry. All arguments will 
be local, except to note implications for conformal infinity, if it exists. 
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II. Geometrical properties 

The line-element may be written locally in double-null form as 

ds 2 = r 2 d0 2 - 2e _/ df +d£- (1) 

where dfi 2 refers to the unit sphere and r and / are functions of the null coordinates 
(£ + , This double-null form is natural in the sense that each symmetric sphere has two 
preferred normal directions, the null directions d/d^. One may also use one spatial and 
one temporal direction, as in Section IV, but there is no unique choice of such directions, 
which makes it more difficult to check coordinate invariance. In double-null form, the 
remaining coordinate freedom consists simply of the diffeomorphisms 

(2) 

The geometry is given by the metric (1) modulo the null rescalings (2). Note that r is a 
geometrical invariant but that / is not. The area of a symmetric sphere is 4-7rr 2 , so that r 
is the areal radius, and will simply be called the radius. One may take r > 0, with r = 
being discussed shortly. 

The space-time will be assumed time-orientable, and d/d^ will be assumed future- 
pointing. The expansions may be defined by 

0± = 2r~ 1 d±r (3) 

where d± denotes the coordinate derivative along £ ± . The expansions measure whether the 
light rays normal to a sphere are diverging (9 > 0) or converging (9 < 0), or equivalently, 
whether the area of the spheres is increasing or decreasing in the null directions. Note that 
the signs of 9± are geometrical invariants, but their actual values are not. An invariant 
combination is e-^6>+6>_, or equivalently g~ 1 (dr, dr) = — |e^r 2 6>+6>_, where g is the space- 
time metric. Indeed, the only invariants of the metric and its first derivatives are functions 
of r and 9+9-. The latter invariant has an important geometrical and physical meaning: 
a metric sphere is said to be (i) trapped if 9+9- > 0, (ii) untrapped if 9 + 9- < and (iii) 
marginal if 9+9- = 0. (Equivalently, if (? _1 (dr) is temporal, spatial or null respectively). 
This terminology for surfaces will be extended to hypersurfaces and space-time regions. If 
9+9- is a function with non-vanishing derivative, the space-time is divided into trapped 
and untrapped regions, separated by marginal hypersurfaces. The following subdivisions 
may be made. 

(i) A trapped sphere is future if 9± < and past if 9± > [6-7]. (Equivalently, if g _1 (dr) 
is future-temporal or past-temporal respectively). Future and past trapped spheres 
occur in black and white holes respectively, and also in cosmological models. 

(ii) On an untrapped sphere, a spatial or null normal vector z is outgoing if z(dr) > 
and ingoing if z(dr) < 0. Equivalently, fixing the orientation locally by 9+ > and 
9- < 0, z is outgoing if g(z, d + ) > or g{z,d-) < and ingoing if g(z, d + ) < or 
g(z,d-) > [8]. In particular, d + and cL are respectively the outgoing and ingoing 
null normal vectors. It is easily checked that the area is increasing in any outgoing 
spatial or null direction, and decreasing in any ingoing spatial or null direction. 
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(iii) A marginal sphere with 9 + = is future if 6>_ < 0, past if 6>_ > 0, bifurcating if 6>_ = 0, 
outer if d-9 + < 0, inner if d-9 + > and degenerate if d-9 + = [6-7]. The closure 
of a hypersurface foliated by future or past, outer or inner marginal spheres is called 
a (non-degenerate) trapping horizon [6-7]. Outer trapping horizons possess various 
easily proven properties which are often intuitively ascribed to black or white holes, 
including confinement of observers and analogues of the zeroth, first and second laws 
of thermodynamics [6-9]. Inner trapping horizons include cosmo logical horizons as 
well as the possible inner boundaries of black and white holes. In this article, future 
(respectively past) outer trapping horizons will be taken as the definition of black 
(respectively white) holes. This enables the discussion of black holes in general space- 
times, not just those which are asymptotically flat. In particular, there will be no 
discussion of event horizons or apparent horizons, which are defined in asymptotically 
flat space-times only [10-11]. The trapping horizon is a more useful concept because it 
is sufficiently general, is defined quasi-locally and has the properties mentioned above. 

The Misner-Sharp spherically symmetric gravitational energy, or simply the energy, may 
be defined in units G = 1 by 

E = \r (1 -^ _1 (dr,dr)) = \r + e* rd+rd-r = \r + \e f r 3 9+9-. (4) 

The form actually given by Misner & Sharp is derived in Section IV. Note that E is an 
invariant. Indeed, the only invariants of the metric and its first derivatives are functions 
of r and E, as explained above. This makes r and E natural variables to use, as has been 
rediscovered many times by different authors. It transpires that remarkably many key 
geometrical properties of spherically symmetric space-times are controlled by r and E. 

Proposition 1: trapping. A metric sphere is trapped if and only if E > |r, marginal if and 
only if E = ^r, and untrapped if and only if E < |r. 

Proof. By definition. □ 

This property is mathematically trivial given the definition in the above form, but is 
physically important because it shows that the ratio E/r controls the formation of black 
and white holes, and trapped spheres generally. Note that the material (or passive) mass 
does not have this property; the sharpest relations [12] between trapped spheres and the 
material mass fall short of necessary and sufficient conditions. In other words, it is not 
the material mass which directly controls the formation of black and white holes, but the 
effective energy E. 

Consider the two-dimensional space-time obtained by taking the quotient by the 
spheres of symmetry. If r = coincides with a boundary of the quotient space-time, then it 
will be called a centre. A central point p will be called regular if g~ 1 (dr, dr) — 1 = 0(r 2 ) as 
p is approached, and singular otherwise. Then E = (9(r 3 ) at a regular centre. A singular 
centre will also be referred to as a central singularity. 

A point in the centre will be said to be trapped if it surrounded by a neighbour- 
hood of trapped spheres, and untrapped if it surrounded by a neighbourhood of untrapped 
spheres. A regular centre is untrapped, but central singularities may be either trapped, 
as in the positive-energy Schwarzschild solution, or untrapped, as in the negative-energy 
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Schwarzschild solution. Whether singularities are trapped or untrapped is relevant to Pen- 
rose's cosmic censorship hypothesis [13-14]. Indeed, one might formulate a local version of 
weak cosmic censorship [13] in terms of whether singularities are trapped. Directly relevant 
to strong cosmic censorship [14] is whether the singularity is causal or spatial, defined with 
respect to the quotient metric. It turns out that both features are controlled by E [15-16]. 

Proposition 2: central singularities. For a central singularity at p, if E > E in a neigh- 
bourhood of p for some constant E > (respectively E < E < 0) then: (i) p is trapped 
(respectively untrapped); (ii) if the singularity is differentiable at p, then it is spatial 
(respectively temporal) . 

Proof. Consider the case E > E > (the case E < E < being similar), (i) In the 
neighbourhood, ^e^r 3 9+9- = 2E — r > 2Eq — r > for sufficiently small r < 2Eq. So 
9+9- > in a neighbourhood of p, i.e. this neighbourhood consists of trapped spheres, (ii) 
The tangent vector z = d/d( to the singularity is a linear combination z = (3d+ — ad- of 
the null normals d±, so that = dr/d( = f3d + r — ad-r. But d + rd-r = \r 2 9 + 9- > in 
the above neighbourhood, so a/3 > 0, which means that z is spatial. □ 

This reflects a physically important idea: there is a connection between the sign of the 
energy, the causal nature of singularities and whether they are trapped. Combined with 
the property E > 0, which will be derived under certain assumptions in Proposition 6, 
this result supports the cosmic censorship hypothesis. It falls short of a proof of cosmic 
censorship for two reasons. Firstly, it is possible to have E < when the assumptions of 
Proposition 6 do not hold, as will be discussed at the end of Section III. Secondly, there 
is the case where E\ p = 0, in which case p could be a spatial, null or temporal singularity, 
or a regular centre. Specifically, if 2E/r\ p > 1 the singularity is spatial and trapped, 
and if 2E/r\ p < 1 the singularity is temporal and untrapped. If 2E/r\ p = 1, one must 
look at higher orders, (2E/r — l)/r, and so on. Exactly which possibility occurs seems to 
depend on the matter field. According to the analysis of Christodoulou [17], for a massless 
scalar field it is possible to obtain causal central singularities, but such configurations are 
non-generic with respect to initial data. Conversely, for pure radiation (or null dust), a 
sufficiently weak wave travelling into an initially flat space-time necessarily creates a null 
singularity which is at least locally visible [18-19]. Such visible null singularities are also 
possible for dust [20-22]. 

Despite such material-dependent differences, one useful fact remains: if E > 0, a 
central singularity which is either causal or untrapped must have vanishing energy, E\ p = 0. 
(For instance, the analysis of Joshi & Dwivedi [23] concerns such singularities). This 
at least constrains counter-examples to cosmic censorship. Moreover, one might expect 
that zero-energy singularities are unstable or non-generic in some sense. Similarly, it has 
been suggested [24] that zero-energy singularities are non-gravitational and do not conflict 
with the spirit of the cosmic censorship hypothesis. Beyond spherical symmetry, there 
is also some evidence for a weakened form of cosmic censorship in which positive-energy 
singularities are censored [25]. 

Proposition 3: asymptotics. In an asymptotically flat space-time, E coincides with the 
Bondi-Sachs (scalar) energy Ebs at null infinity, and with the Arnowitt-Deser-Misner 
(scalar) energy Eadm at spatial infinity. 
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Proof. By definition [2-3,11,15,26,27]: 




(5) 



where 1Z is the Ricci scalar, /j, the area form and A = J ji the area of a family of affinely 
parametrised surfaces lying in a null hypersurface approaching . In spherical symmetry, 
/ n = 4nr 2 and TZ = 2/r 2 , so that 



Ebs = nm E 



(6) 



7^P 



for a null curve 7 approaching p e T^. Similarly, by definition [1,15,28]: 



= nm — 




(7) 




surface approaching z , where a± are the shears corresponding to the expansions 9±. In 
spherical symmetry, a± = 0, so that 



for a spatial curve 7 approaching p = i°. □ 

The result shows that the asymptotic energies (in spherical symmetry) are just special 
cases of E, defined at infinity in an asymptotically flat space-time. It is usual to interpret 
the asymptotic energies as measuring the total energy (or mass — see Appendix C) of an 
asymptotically flat space-time, whereas E provides a more general definition of energy 
which applies locally as well as asymptotically. 

III. Dynamical properties 

Having derived various purely geometrical properties of E, consider now applying the 
Einstein equations. The most general form of the Einstein tensor in spherical symmetry 
determines the most general material stress-energy tensor T, given by the line-element 



E adm = lim E 



(8) 



7^P 



pr 2 d0 2 + T ++ (d£+) 2 + T__(dC ) 2 + 2T + _d£ + d£-. 



(9) 



The Einstein equations are 



d±d±r + d±fd±r = —4i:rT±± 
rd+d-r + d+rd-r + \e~ } = 4nr 2 T + _ 

r 2 d + d-f + 2d + rd-r + e~ f = 8yrr 2 (T + _ + e~ f p) 



(10a) 
(106) 
(10c) 
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recalling the units G = 1. The variation of E is determined by these equations as 

d±E = 4ne f r 2 (T+_d±r - T±±d T r) = 2ne f r 3 (T + _9± - T±±9 T ). (11) 
This can also be written in a manifestly covariant form [29-30] . 

Proposition 4-' vacuum. In vacuo, E is constant and the solution is locally isometric to the 
Schwarzschild solution with energy E. 

Proof. In vacuo, d±E = 0, so E is constant. The vacuum Einstein equations are 

d±(e f d ± r)=0 (12a) 
d+d-r = -Ee- f r~ 2 (126) 
d+d-f = -2Ee- f r~ 3 . (12c) 

If E ^ 0, a straightforward calculation yields <9_|_<9_(/ — logr — r/2E) = 0, with general 
solution / = log r+r/2E + f + + f~ for some functions / ± (^ ± ) of integration. The rescaling 
freedom (2) corresponds to choice of / ± , so that one may fix coordinates such that 

e f = re r/2E /16E 3_ (13) 

The equations (12a) integrate to e?d±r = rf^ for some functions ?y ± (£ ± ) of integration. 
Differentiating using (12b), d±i] ± = —1/4E, which integrates to ^ = —t^/AE, fixing the 
zero. Thus 

= (1 - r /2E)e r / 2E (14) 

which implicitly determines r (and hence /) as a function of . This is the Kruskal 
form of the Schwarzschild solution with energy E, which can be put in static form in terms 
of t^ 1 = 2E\og(—^ ± /^ Zf ). Similarly, if E = then flat space-time is recovered. □ 

This is a proof of Birkhoff 's theorem: a vacuum, spherically symmetric space-time must 
be the Schwarzschild solution. The proof improves on the usual one [10] in that a global 
coordinate chart is obtained automatically, so that one does not have to subsequently join 
the r > 2E and r < 2E Schwarzschild charts. 

The Schwarzschild solution provides an example of Propositions 1-3. For the case 
E > there are trapped spatial singularities, while for the case E < there is an untrapped 
temporal singularity. In the former case, there are trapped spheres in the black-hole and 
white-hole regions E > |r, and the event horizons coincide with the trapping horizons 
E=\r. 

Consider now the non- vacuum cases. In order to obtain results which are as general as 
possible, the type of matter will not be fixed but energy conditions will be imposed instead. 
Three useful energy conditions are as follows [10-11]. The null energy (or convergence) 
condition states that a "null observer" measures non- negative energy: 

NEC: g{u,u) = T(u,u)>0. (15) 

The weak energy condition states that a causal observer measures non-negative energy: 

WEC: g(u,u)<0 T(u,u)>0. (16) 
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The dominant energy condition states that a future-causal observer measures future-causal 
momentum: 



DEC: 



g(u, u) < 0, g(v, v) < 0, g(u, v) < T(u, v) > 0. 



(17) 



Clearly 



DEC WEC NEC. 



(18) 



All this applies to general space-times. In the spherically symmetric case, 



NEC =>■ T±± > 



(19) 



and 



DEC =>■ T±± > 0, T + _ > 0. 



(20) 



Proposition 5: monotonicity. If the dominant energy condition holds on an untrapped 
sphere, E is non-decreasing (respectively non-increasing) in any outgoing (respectively 
ingoing) spatial or null direction. 

Proof. Fix the orientation of the untrapped sphere by 9 + > and 6>_ < 0. The variation 
formula (11) and dominant energy condition (20) yield d + E > and d-E < 0, i.e. E 
is non-decreasing (respectively non-increasing) in the outgoing (respectively ingoing) null 
direction. If z = d/d( is an outgoing spatial vector, then z = (3d+ — ad- with a > 
and j3 > 0, which yields dE/d( = (3d+E — ad-E > 0. Similarly for an ingoing spatial 
direction. □ 

The proof illustrates the economy of the double-null approach: monotonicity in any spatial 
direction follows immediately from monotonicity in the null directions. This monotonicity 
property has the physical interpretation that the energy contained in a sphere is non- 
decreasing as the sphere is perturbed outwards. Note that the result is for untrapped 
spheres only, though a similar result for marginal spheres will be given in Proposition 8. 
There is no possibility of a similarly general monotonicity result for trapped spheres, since 
they do not have a locally preferred orientation. 

Proposition 5 A: asymptotic monotonicity. If the dominant energy condition holds at X + 
(respectively X~), then Ebs is non-increasing (respectively non-decreasing) to the future. 

Proof. By Propositions 3 and 5, 8-Ebs = <9_(lim 7 ^ p E) = \\va 1 ^ p d-E < 0. Here 
commutativity follows from the asymptotic expansions [11], as in [8]. □ 

At X + , this is the Bondi-Sachs energy-loss property, which is usually interpreted as describ- 
ing a loss of energy due to outgoing radiation. Similarly, the more general monotonicity 
property of E may be interpreted as being due to ingoing and outgoing radiation. 

Positivity of E is already known in some circumstances [30] but can be established 
directly and quite generally from the monotonicity property. Note first that E may be 
negative, as in the negative-energy Schwarzschild solution. On the other hand, Proposition 
1 shows that E is automatically positive for trapped and marginal spheres, so that it is 
only for untrapped spheres that E might be negative. Positivity on untrapped spheres 
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is established below in two physically relevant contexts, namely where there is a regular 
centre and where there is a black or white hole. 

Proposition 6: positivity. If the dominant energy condition holds on an untrapped spatial 
hypersurface with regular centre, then E > on the hypersurface. 

Proof. By Proposition 5, since E = at a regular centre. □ 

Moreover, denoting the tangent vector to the hypersurface as before by d/d( = (3d + — ad-, 
with C = being the centre, one can write explicitly 

E(C) =4tt [ C e f r 2 {(f3T+_ + aT__)d+r - (aT+_ + f3T ++ )d-r}d(. (21) 
Jo 

The positivity property has the physical interpretation that under the stated circum- 
stances, total energy cannot be negative. This is not immediately obvious even given an 
energy condition on the matter, since gravitational potential energy tends to be negative. 
The result shows that the total energy E, including potential energy, cannot be negative. 

Proposition 7: area inequality. If the dominant energy condition holds on an untrapped 
spatial hypersurface bounded at the inward end by a marginal sphere of radius ro, then 
i? > |r on the hypersurface. 

Proof. By Propositions 1 and 5. □ 

Since r > 0, this is a stronger result than mere positivity of E: there is a positive lower 
bound on E. The physical interpretation is that if there is a black or white hole of area 
4-7rrQ, then the energy measured outside the hole is at least ^ro- As for the positivity 
result, one can write an explicit formula 

E(C) = \r + 4tt [ C e f r 2 {(f3T+_ + aT__)d+r - («T+_ + /3T ++ )d_r}d( (22) 

J Co 

where C = Co is the marginal surface. 

Proposition 6 A: asymptotic positivity. If the dominant energy condition holds on a spatial 
hypersurface which has a regular centre and extends to (respectively then Ebs > 
(respectively Eadm > 0) there. 

Proof. If the hypersurface is untrapped, Propositions 3 and 6 suffice. Otherwise, the 
hypersurface contains an outermost marginal sphere, so Propositions 3 and 7 suffice. □ 

This combines the famous positive-energy theorems for the spherically symmetric case. 
Note that the energy can be negative if the above assumptions do not hold. For instance, 
the asymptotic energies are negative for the negative-energy Schwarzschild solution. 

Proposition 7 A: asymptotic area inequality. If the dominant energy condition holds on a 
spatial hypersurface which contains an outermost marginal sphere of radius ro and which 
extends to J ± (respectively then E B s > \fQ (respectively Eadm > \tq) there. 

Proof. By Propositions 3 and 7. □ 
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The result is the spherically symmetric case of the isoperimetric inequality conjectured 
by Penrose [5], who argued that it was required by the cosmic censorship hypothesis. 
Establishing the inequality even in spherical symmetry appears to be new. It was recently 
established for maximal hypersurfaces in spherical symmetry [31]. 

The properties 5A, 6A and 7A of the asymptotic energies are of interest in their own 
right. Nevertheless, they are just special cases of properties of E. If these properties 
of the asymptotic energies are accorded their usual conceptual and physical importance, 
then the more general properties of E are of even greater importance. The idealisation of 
asymptotic flatness is no longer necessary for the formulation of such ideas about energy. 

Proposition 8: second law [6]. If the null energy condition holds on a future (respectively 
past) outer trapping horizon, or on a past (respectively future) inner trapping horizon, 
then E = is non-decreasing (respectively non-increasing) along the horizon. 

Proof. Denote the tangent to the horizon by d/d( = f3d + — ad- and fix the orientations 
by 6+ = and /3 > on the horizon. Then = d6+/d( = f3d + 6 + - ad-0 + yields 
dr/d( = —ad-r = —f3r6-d + 6 + /2d-6 + . The focussing equation (10a) and null energy 
condition (19) yield d+9+ < 0, and the signs of 6>_ and d-9+ are given by the definition of 
future or past, outer or inner trapping horizons. Thus dr/d( > for future outer or past 
inner trapping horizons, and dr/d( < for past outer or future inner trapping horizons. 
□ 

Propositions 5-8 may be loosely summarised as follows. If the dominant energy condi- 
tion holds, E is non-decreasing in outgoing directions, defined for untrapped or marginal 
surfaces, including at conformal infinity if it exists. So in the untrapped region outside 
a regular centre, E is non-negative, including at conformal infinity if it exists for the un- 
trapped region. But E is also positive in a trapped region. Moreover, in the untrapped 
region outside a black or white hole (defined by outer marginal surfaces) of area 47rrg, 
E > ^ro, including at conformal infinity if it exists for the untrapped region. 

Propositions 1-8 alone give a quite coherent picture of gravitational collapse, which 
may be further refined using related results [6-9] . Recall first the expected picture accord- 
ing to the cosmic censorship hypothesis [13-14]. Here a centre which is initially regular 
subsequently becomes singular, but in such a way that the singularity is generically spa- 
tial and trapped. That is, a future trapping horizon forms before the singularity. More 
precisely, one expects the trapping horizon to intersect the centre, typically at the first 
singular point. f The horizon develops outwards with non-decreasing E = |r. Outside the 
horizon is an untrapped region in which < E < \r. Inside the horizon is a trapped region 
in which E > |r, which is expected to extend to the centre r = in such a way that E is 
generically positive there. This picture has been confirmed for the massless scalar field by 
Christodoulou [17]. 



f This does not contradict the fact that, when using a family of spatial hypersurfaces 
developing in time, trapped surfaces are often first detected away from the centre. This 
just depends on the choice of hypersurfaces; if the horizon is spatial, one may choose a 
spatial hypersurface to touch it anywhere. There is no unique first moment of trapped 
surface formation, except in a limiting sense at the centre. 
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To what extent do the above results support this picture? Firstly, there is the connec- 
tion of Proposition 2 between the causal nature of central singularities and whether they 
are trapped. Thus two aspects of cosmic censorship are linked. Secondly, both aspects are 
linked to the sign of the energy E, again a physically important connection. Thirdly, there 
are the positivity properties of E, Propositions 1, 6 and 7. Consider then a space-time with 
partially regular centre, and the external region X of points connectable to the regular 
part of the centre by a non-temporal curve (respecting the spherical symmetry). As long 
as X remains untrapped, E is positive and no central singularities can form. If part of 
X is trapped, i.e. a trapping horizon forms, then inside the trapped region any central 
singularity must be either spatial or have zero energy. One expects that zero-energy singu- 
larities are non-generic or unstable in some sense. The only other way that a non-spatial 
central singularity can form in X is if a second trapping horizon forms, separating the 
trapped region from another untrapped region. That is, the black or white hole must have 
an inner boundary. This is also possible, as the Reissner-Nordstrom solution shows, but 
again such horizons are thought to be unstable by the energy-inflation effect [32-34]: for 
a test field on such a background, the perturbation in E typically becomes unbounded 
at the horizon. This leaves the only other possibility for non-spatial central singularities 
as being outside X. Again this is possible: the centre may simply become singular but 
remain causal, implying that E is no longer positive. There are such examples [17-24], 
but it is noteworthy that they have zero rather than negative E. This suggests that there 
may be some mechanism forbidding such formation of negative-energy singularities. As to 
the zero-energy singularities, one might again expect them to be unstable in some sense. 
Both possibilities indicate the need to study the behaviour of E as the centre becomes 
singular. Finally, there is the possibility of non-central singularities, such as shell-crossing 
singularities [35-36]. However, these have been found only for matter fields which admit 
similar singularities as test fields on flat space-time, and therefore may be dismissed as 
pathological matter models. 

To conclude: to find a convincing counter-example to cosmic censorship, i.e. the stable 
formation of a causal singularity, one needs to study the possibility of negative-energy 
singularities, the stability of zero-energy singularities or the stability of inner horizons. It 
is noteworthy that E plays a crucial role in each case. The importance of energy for cosmic 
censorship seems clear. 

IV. Spatial hypersurfaces and the Newtonian limit 

Consider any spatial hypersurface E. Set up coordinates (r, Q such that d/d( is tangent 
to E and d/dr is orthogonal to d/d( : with r being proper time: 

g(d/dr,d/d() = (23a) 
g(d/dr,d/dr)= -1. (236) 

Define a function A by 

e x =g(d/dCd/d(). (24) 
For any such coordinates (r, £), the rescaling freedom (2) can be used to fix ^ such that 

v / 2e-^ 2 dC ± = dr ± e A / 2 dC (25) 
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The line-element (1) transforms to 



ds 2 = r 2 d0 2 + e A dC 2 - dr 2 (26) 

where r and A are functions of (r, £). Denote <^> = d(p/dr and <// = d(p/d(. Then the 
definition (4) of E can be rewritten as 

1- — = e -V) 2 -r 2 (27) 
r 

which is the form actually given by Misner & Sharp [4], with different notation. 

The Einstein equations may be transformed to these coordinates, but for the following 
results it suffices to find the corresponding variation formulas for E, which are 

E' = 47rr 2 (T 00 r' — T 01 r) (28a) 

E = 47rr 2 e" A (T ir' - Tur) (286) 

where 

Too = T(d/dr, d/dr) (29a) 

T 01 = T(d/dr,d/d() (296) 

T 11 =T(d/dCd/d(). (29c) 

Proposition 9: small spheres. Near a regular centre with tangent d/dr, if T o = O(l) and 
T i = 0(1) then 

E= |7rr 3 T 00 | r =o + 0(r 4 ). (30) 

Proof. Near a regular centre r = 0, r' = 0(1). If d/dr is tangent to the centre, then 
r = O(r). Thus E' = 4irr 2 Toor' + 0(r 3 ), which integrates along the hypersurface E to the 
above result. □ 

In other words, the leading term in E is the product of volume |7rr 3 and density T o, as 
would be expected physically. 

Misner & Sharp [4] derived a useful formula for E which will be derived below and 
used to find the Newtonian limit. Misner & Sharp considered a perfect fluid with energy 
density p, pressure p and velocity d/dr: 

T = (p + p)dr ®df + pg (31a) 
g(d/dr,d/dr) = -l. (316) 

Adapt the (r, Q coordinates to the fluid by taking r = f. Then the variation formulas 
(28) for E reduce to 

E' = 4nr 2 pr (32a) 

E = -4nr 2 pf. (326) 
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The second equation (32b) expresses the rate E of work done by the force 47vr 2 p due to 
pressure. In the absence of pressure, E is conserved, E = 0. The first equation (32a) 
allows E to be expressed as an integral of e over E. Specifically, if E has a regular centre 
C = 0, (32a) integrates to 

E(Q= f 47rrV'dC- (33) 



o 



Since the volume form *1 of E is given by (26) as 



*ip - 
t Jo 



47rr 2 e A/ VdC (34) 



this may be rewritten as 



E 



1 /2 

= J *pe~ x / 2 r' = J *p(\ + r 2 -^ (35) 



where the second expression follows from (27). Comparing with the material (or passive) 
mass 



M = *T 00 = *P (36) 
Jt Jt 

it can be seen that the integrand for E differs from that of M by a factor which Misner & 
Sharp interpreted as being due to kinetic and potential energy. This can be made precise 
in the Newtonian limit in terms of the kinetic energy K and gravitational potential energy 
V, defined as explained in Appendix A by 

K= [ *\pf 2 (37a) 
V=-f*¥*. (376) 

Factors of the speed of light c may be introduced on dimensional grounds by the formal re- 
placements r h-> cr, (r, *1) h-> (r, *1), (p, M) i-> c~ 2 (p, M), (p, V, E 1 ) i-> c~ 4 (p, V, E 1 ), 
assumed in the remainder of this Section. These factors are determined simply by the 
desired interpretation of the various quantities as time, length, mass etc. 

Proposition 10: Newtonian limit. For a perfect fluid on a spatial hypersurface E with 
regular centre: if (*1, r, f, p) = 0(1) as c — > oo, then (M, if, V) = 0(1) and 

E = Mc 2 + if + y + 0(c- 2 ). (38) 

Additionally, if p = 0(1) then Newtonian conservation of mass is recovered: 

M = 0( C - 2 ). (39) 
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Proof. Inserting the factors of c, (36) and (37) take the same form, so that (M, K, V) = 
0(1), while (35) takes the form 




Thus E = c 2 J E *p + 0(1) = Mc 2 + 0(1) to leading order. Expanding the square root, 
E = Mc 2 + / s *p(\f 2 - E/rc 2 ) + 0(c" 2 ) = Mc 2 + K + V + 0(c" 2 ). Similarly, (32b) 
expands to -4%r 2 pf = Mc 2 + 0(1). □ 

In words, E yields the Newtonian mass M to leading order and the Newtonian kinetic 
energy K and gravitational potential energy V to next highest order. This illustrates 
how E measures the combined energy including contributions from mass, kinetic energy 
and potential energy. Note also that the quantities M , K and V are all defined in the 
full theory rather than just the Newtonian limit. Of these, M may always be interpreted 
as the material (or passive) mass, whereas the interpretation of K and V as kinetic and 
potential energy is reliable in the Newtonian limit only. In general, E cannot be expressed 
as a sum of individually meaningful energies, as the form (35) indicates. It is only the 
combined energy E which is physically meaningful in general. 

Incidentally, the above considerations provide one reason for referring to E as an 
energy rather than a mass, though the latter is more common. Although mass and energy 
are formally equivalent in Relativity, the two words carry connotations inherited from their 
status as distinct concepts in Newtonian theory. Specifically, mass is a measure of matter 
whereas energy exists in many different forms. So it is reasonable to describe M as mass, 
since it is simply the integral (36) of the material density, and to describe E as energy, 
since it contains contributions from kinetic and potential energy. 

Under the assumptions of Proposition 10, the inverse metric is Euclidean to leading 
order, g~ x = 5 _1 + 0(c -2 ), where the flat metric 6 is given by r 2 d0 2 + dr 2 , so that r is a 
standard radial coordinate. Thus the flat space of Newtonian theory is recovered. 

Note that it has not been necessary to introduce the usual coordinate conditions [37] 
required to obtain an inertial frame in which Newton's laws hold, or to obtain Newtonian 
solutions from solutions to the Einstein equations [38]. In contrast, the Newtonian be- 
haviour of E may be obtained simply by expanding a formula for E in powers of c. In 
this sense, the recovery of Newtonian mass and energy from E is more robust than the 
complete recovery of Newton's gravitational theory from Einstein's theory. This could be 
interpreted as meaning that energy is a fundamental concept which connects Newtonian 
theory with Relativity. Certainly this is consistent with the key role that the equivalence 
of mass and energy played in the historical development of Relativity. 

V. Test particles and the special-relativistic limit 

The energy of a spherical shell of test particles may be defined as the variation in E 
induced by the perturbation of the test particles, as suggested by Hajicek [39]. Denote 
the perturbations of (T,E,M) by (T, £,A4). In terms of the constant rest mass m and 
velocity d/df of the shell, the stress-energy tensor is given analogously to (31) by 

T = m5dr <g> df (41) 
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where 5 is the Dirac distribution with support at the intersection S of the shell with a 
spatial hyper surface E: 

J^*5cp = ip\ s . (42) 

Evaluation on S is clear in context and will be omitted henceforth. Take E to be orthogonal 
to the shell and adopt the notation of Section IV for spatial hypersurfaces, identifying 
r = f. Then the Hajfcek energy £ is defined by the perturbation of (11) or (28), 

£' = 4nr 2 (T 00 r' - T 01 r) (43a) 

£ = 47rrV A (T 01 r' -T n r). (436) 

and M by (36), 

M = f *T 00 . (44) 

Proposition 11: test particles. For a spherical shell of test particles, (i) the perturbations 
in M and E are respectively 

M = m (45a) 



8 = ±m [ 1 + r 2 - — ) (456) 



2E 



1/2 



where the sign is that of r'. Both are conserved: 



M = (46a) 
8 = 0. (466) 



(ii) In the Newtonian limit, with assumptions as for Proposition 10, 



±8 = mc 2 + \mf 2 - — + 0(c~ 2 ). (47) 



Proof, (i) Since T 00 = m5 and T 01 = Tn = 0, (44) reads M = J E *m5 = m, (43b) 
reads 8 = and (43a) reads £' = A-nr 2 r'm5, which integrates using (27) and (34) to 
£ = ± / E *mS(l + r 2 - 2E/r) 1 / 2 = ±m(l + r 2 - 2^/r) 1 / 2 . (ii) Inserting the factors of c, 

= + V2 < 48 > 

which expands as above. □ 

In words: the perturbation in M is just the rest mass m, while the perturbation in E 
contains contributions from energy, in particular the correct Newtonian expressions for 
kinetic energy \mf 2 and gravitational potential energy —mM/r. Combined with the fact 
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that £ is conserved as the shell evolves, this confirms the physical interpretation of £ as 
the active energy of the shell of test particles. 

Conservation (46b) of S, written explicitly using (45b), is equivalent to the equation 
of motion of the test particles: 

r = -$/r (49) 

where $ = —E/r. This has a similar form to the Newtonian equation of motion of test 
particles in a central gravitational field, with E playing the role of the active gravitational 
energy of the source and $ the gravitational potential. Again this confirms the physical 
interpretation of E as the active gravitational energy. 

Proposition 12: special-relativistic limit. In flat space-time, the perturbation in E for a 
spherical shell of test particles is 

S\ E=0 = ±m(l + r 2 ) 1/2 . (50) 

Proof. By Proposition 11. □ 

This has been dignified as a separate Proposition because it is important that E yields 
the correct result in the special-relativistic limit, namely that £ reduces to the standard 
expression (50) for the energy of test particles in special relativity. In particular, £ is not 
just the rest mass. This is noteworthy since the opposite is sometimes claimed, namely 
that E corresponds to the rest mass rather than energy of a test particle in flat space- 
time. This misunderstanding seems to be the reason — other than apathy — that E is usually 
described as mass rather than energy. On the contrary, using the terminology standard 
in both special relativity and Newtonian theory, £ and E are energies and m and M 
masses. In the same sense, the following are all energies rather than masses: quasi-local 
energies which agree with E in spherical symmetry, such as those of Penrose and Hawking 
[15,40-46]; the asymptotic (scalar) energies Eadm and Ebs', the Schwarzschild energy. 

The mistaken claim that E corresponds to m rather than £ is sometimes ascribed to 
the fact that, in flat space-time, £ depends on the choice of frame (or observer) whereas m 
is frame-independent, like E. But except in flat space-time, £ is also frame-independent: 
the velocity f is relative to the preferred frame determined by r, not the frame of some 
observer. Only in flat space-time is r not unique, due to the Lorentz symmetry. It is 
this degeneracy of flat space-time that renders £ frame-dependent in that special case. 
Frame-dependence is not a general feature of energy, as discussed further in Appendix C. 
For instance, Newtonian gravitational potential energy is frame-independent. 

VI. Remarks: gravitational energy in general 

The Misner-Sharp energy E has an impressive variety of useful properties, ranging from the 
asymptotically flat, vacuum, small-sphere, Newtonian, test-particle and special-relativistic 
limits to the black holes and singularities characteristic of strong gravitational fields. These 
properties have an exact geometrical character and are simultaneously of direct physical 
relevance. In particular, E has quite general monotonicity and positivity properties, de- 
termines the causal nature of central singularities, and determines when trapping occurs. 
This makes E useful in many different applications, particularly regarding black holes and 
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singularities. Indeed, E has been rediscovered by many authors and is often crucial to 
their analyses [16-22,24,29-30,32-34,36,39]. In particular, E arises as the charge associ- 
ated with a conserved current, as explained in Appendix B, and can be derived from the 
"2+2" Hamiltonian [15]. One might also speculate that the radius r and energy E are rel- 
evant in quantum theory [47] . Certainly E plays a key role in the vacuum case, according 
to the analysis of Kuchaf [48] . 

The conceptual and physical importance of E in spherical symmetry encourages the 
search for a more general definition of gravitational energy. It is widely accepted that such 
a general definition should reduce to E in spherical symmetry, as occurs for definitions in 
the style of Penrose [40-44] or Hawking [15,45-46]. Exceptions include the Brown- York 
energy [49], which gives a value different from E for the Schwarzschild solution, and the 
Bartnik energy [50] , which is undefined for trapped surfaces. Whether such definitions have 
some other physical meaning is unclear, but they do not represent gravitational energy in 
the sense of Propositions 1-12. 

Familiarity with the spherically symmetric case also yields other guidelines to more 
general definitions. Firstly, it is noteworthy that E is defined on spheres rather than 
hyper surfaces. One can write an implicit expression for E involving an integral (35) over 
a hypersurface with regular centre, if such a hypersurface exists, but the general definition 
(4) or (27) is a function of spheres only. This is quite different to the definition of mass 
M as an integral (36) over a hypersurface. So rather than looking for a definition of 
gravitational energy as a hypersurface integral, as would be natural in Newtonian theory, 
one should look for a surface integral. Specifically, one wants an invariant of the intrinsic 
and extrinsic curvature of an embedded surface, which following Penrose [40] may be called 
quasi- local energy. 

Another guideline to more general definitions is that one would like analogous pos- 
itivity and monotonicity properties, and analogous relations to black holes and possibly 
to singularities which are central in some sense. The Hawking energy [45] takes the value 
167T on a spherical marginal surface of area A, and so automatically generalises Propo- 
sitions 1 and 8. The Hawking energy also has the same small-sphere behaviour as in 
Proposition 9 [51]. Generalisations of the positivity and monotonicity Propositions 5-7 
for the Hawking energy have also been found [8]. A warning should be sounded that such 
positivity and monotonicity theorems involve certain assumptions without which the con- 
clusions are invalid. Even the Schwarzschild energy may be negative. Thus there is little 
point in searching for a definition of gravitational energy which is always non-negative by 
definition, as is sometimes suggested [50]. 

The final lesson of the spherically symmetric case is that a general definition of grav- 
itational energy should give the correct results in various limits: the asymptotically-flat, 
vacuum, small-sphere, Newtonian, test-particle and special-relativistic limits. In partic- 
ular, one would like to obtain the Newtonian mass to first order in c, with corrections 
interpretable as Newtonian energies. This is possible for the Penrose energy in certain 
cases [52] . This brings the discussion back to the physical meaning of gravitational energy 
as expressed in the Introduction: an effective energy which is measurable on a surface and 
which is produced by the non-local, non-linear interaction of the mass of sources with the 
energy of the consequent gravitational field. 
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Appendix A: Newtonian gravitational energy 

Recall the dynamical role of the kinetic energy K and potential energy V (37) in Newtonian 
gravity, taking the case of dust, p = 0. The Lagrangian L = K — V determines the Euler- 
Lagrange equation 

d fdL\ dL f /.. M\ , . 



dr \ dr J dr Jy \ r 2 J 

regarding M (or equivalently *p) as an independent variable which is conserved, M = 0. 
This yields Newton's inverse-square law of gravitation:| 

The Hamiltonian H = K + V is therefore conserved: 

H= f *pf ( r + % ) = 0. (A3) 



It is remarkable that both conserved quantities with a clear physical meaning in Newtonian 
gravity, namely the mass M and the Hamiltonian energy H, are encoded in E via the 
expansion E = Mc 2 + H + 0(c~ 2 ) of Proposition 10. 

There are two other expressions which are each sometimes suggested as the potential 
energy of the Newtonian gravitational field, namely 

F = — - / *|V$| 2 (A4a) 
8tt Jy. 

W = \ [ *$p (A4b) 
Jy 

where $ is the Newtonian potential. However, neither F nor W coincide with the potential 
energy V except at infinity, as follows. Integrating V by parts yields 



1 f M 2 



dY 87r Jy r 4 



2r 

Similarly, rewriting W using the Poisson equation 

V 2 $ = 4np (A6) 



| An inverse-square law also exists in the Einstein theory for dust, with M replaced 
by E in (A2). This reinforces the interpretation of E as an active energy. In this case, 
M = E = 0, so that Lagrangian and Hamiltonian formulations are possible with M 
replaced by E in the definition of V. 
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and integrating by parts yields 



W = ir 2 $V n $| 9E - -L j *|V$| 2 (A7) 

where n is the outgoing unit vector. The Poisson equation integrates to V$ = nM/r 2 , so 
that 

M 2 

V+— =F = W-iM$| as . (A8) 

9E 



2r 



So V, F and W all differ by boundary terms. In particular, at finite radius, V and F differ 
unless M vanishes, while V and W differ unless M is constant. However, if the total mass 
lim^oo M is finite and the arbitrary constant in $ is fixed as usual by lim^oo $ = 0, 
then all three quantities coincide at infinity: 

lim V = lim F = lim W. (A9) 

r^oo r^oo r^oo 

Thus all three expressions give the correct asymptotic energy. Nevertheless, at a finite 
radius, only V gives the correct potential energy in the Lagrangian or Hamiltonian sense. 
Thus it is V rather than F or W that should occur in the Newtonian limit of a quasi-local 
energy, even though each would suffice for the asymptotic energy. 

Appendix B: conserved currents and charges 

Kodama [53] showed that spherically symmetric space-times have a conserved current and 
defined an energy as the associated charge. This energy is the Misner-Sharp energy, as 
follows. Kodama introduced vectors k and j which may be defined, up to orientation, by 

/c = curlr (SI) 

and 

g(j) = ~T(k) (B2) 
where the curl refers to the two-dimensional quotient space-time. It follows that 

2E 

g(k,k)= — -l (S3) 

so that k reduces in vacuo (Schwarzschild) to the stationary Killing vector. One may 
interpret k as giving a preferred time for the preferred class of constant-radius observers, 
and j as the corresponding energy-momentum density. Note that k is spatial in trapped 
regions and causal otherwise. If the dominant energy condition holds, j is also causal (or 
zero) in untrapped regions. Fixing the orientation so that k and j are future-causal for 
untrapped spheres, whose orientation is chosen as in the main text, one finds explicitly 
that 



k = e / (d+rd- - d_rd+) = e" A/2 (r' ^ - f-j^J (B4a) 
i8 + Ed_- a _E8 + ) = T -^{Ei--E-y (Bib) 
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Thus k is tangent to the constant-r hypersurfaces and j is tangent to the constant-i? 
hyper surfaces: 

k(r) = (B5a) 

j(E) = 0. (Bbb) 

One also has 

4yrr 2 j = cmlE. (B6) 

Since divcurl = and the four-dimensional divergence is given by Divu = r _2 div(r 2 -u), it 
follows from (Bl) and (B6) that both k and j are covariantly conserved: 

DivA; = (B7a) 
Div j = 0. {Bib) 

By Gauss' law, any such conserved current u has a charge Q u = — J E *g(u, v) which is 
independent of E (with fixed boundaries). Here the notation is that of Section IV, with 
unit normal vector v = d/dr. Using (34) and (B4), the charges associated with k and j 
are the areal volume and energy respectively: 

Qk = - [ *g(k,v)= [ C 4nr 2 r'd( = fyrr 3 (B8a) 



o 



Qj = - f *9(J, v)= f *T(k, v) = [ C E'dC = E. (B8b) 

</£ JE JO 

The charge Qj is the definition of energy of Kodama [53]. Regarding k as generating 
preferred diffeomorphisms, one may interpret j as the corresponding Noether current and 
E as the Noether charge [54] . 

Appendix C: energy-momentum 

In an asymptotically flat space-time, there are definitions not just of total energy, but of 
total energy-momentum [1,11,26,27], a vector (or covector) at infinity. Problems with this 
interpretation as energy-momentum are described below. Since this interpretation is often 
explained by analogy with special relativity, that case is described first. 

In special relativity, a test particle with mass m and unit velocity vector u has energy- 
momentum p = mu. In the frame of an observer with unit velocity v, the energy is —g(p, v) 
and the momentum is p + g(p, v)v. Similarly, a test field with density p and velocity u has 
energy-momentum density j = pu, which to an observer has energy density —g(j,v) and 
momentum density j + g(j,v)v. On a flat spatial hypersurface E with unit normal v, the 
total energy E, momentum IT and energy-momentum J are 

E=-J*g(j,v) (Cla) 

n = / *(j + g(j,v)v) (Cib) 

J = / *7 = Ev + IT (Clc) 



E 
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and the total mass is \f—g(J, J). 

In general relativity, vectors are not integrable, so IT and J are ill defined. In contrast, 
E is well denned for given j, but generally depends on the hypersurface E. However, it 
will be independent of E (with fixed boundaries) if j is covariantly conserved, by Gauss' 
law. As explained in Appendix B, spherical symmetry admits such an energy- momentum 
density: the Kodama current j, a preferred contraction (B2) of the stress-energy tensor. 
The corresponding energy E agrees (B8b) with the Misner- Sharp energy. Again, this 
identifies E as energy rather than mass. The total momentum IT and energy-momentum J 
remain ill defined. (One could make sense of them by taking the unit spatial vector outside 
the integral, but the result would depend on E). So proceeding by analogy with special 
relativity yields, independent of hypersurface, a total energy but no total momentum, 
energy-momentum or mass. 

Another useful asymptotic quantity is 

K = lim k (C2) 

r — »oo 

where k is the Kodama vector (Bl). By (B3), K is a unit vector, so it may be interpreted 
as the preferred asymptotic velocity, or the total velocity of the source. This provides a 
preferred frame at infinity, the asymptotically stationary frame. 

Now suppose that one believed that E was mass rather than energy. Then by analogy 
with test particles in special relativity, one might define 

P = lim Ek (C3) 

and call it asymptotic energy-momentum. This applies to both null and spatial infinity, 
agreeing with the Bondi-Sachs [11,26,27] and Arnowitt-Deser-Misner [1] energy-momentum 
respectively: P = EbsK and P = EadmK respectively, by Proposition 3. Pursuing this 
idea, one might refer to —g(P, v) and P + g(P, v)v as energy and momentum respectively, 
and to E B s or Eadm as mass. This terminology will be avoided here because E is energy 
rather than mass, according to the Newtonian, test-particle and special-relativistic limits 
of Sections IV and V, as well as the above analogy with test fields in special relativity. (If 
necessary for clarity, Ess or Eadm may be referred to as the scalar asymptotic energy). 
By the same reasoning, P is the product of velocity and energy rather than velocity and 
mass, so is not energy-momentum. 

To measure the gravitational field experimentally, one may observe the motion of test 
particles. The equation of motion (49) shows that the measurable energy is just E, not 
something that agrees asymptotically with —g(P, v). 

If P did measure asymptotic energy-momentum, one might expect it to be the limit of 
a quasi-local energy-momentum, as for the energy E itself. But Ek is not recognisable as 
an energy-momentum at finite radius. Independently, the energy E and current k are both 
useful quantities with simple physical interpretations, but their product has no obvious 
physical meaning. Asymptotically one has the limits lim^oo E and lim^oo k, but their 
product P also has no obvious physical meaning. So it seems that there is no asymptotic 
(or quasi-local) energy-momentum. Instead one has the radius r, the current k, the local 
energy-momentum density j and the quasi-local energy E. 
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